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ABSTRACT 



The discovery of transiting planets around bright stars holds the potential 
to greatly enhance our understanding of planetary atmospheres. In this work 
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we present the search for transits of HD168443b, a massive planet orbiting the 
bright star HD 168443 {V = 6.92) with a period of 58.11 days. The high eccen- 
tricity of the planetary orbit (e = 0.53) significantly enhances the a-priori tran- 
sit probability beyond that expected for a circular orbit, making HD 168443 a 
candidate for our ongoing Transit Ephemeris Refinement and Monitoring Survey 
(TERMS). Using additional radial velocities from Keck-HIRES, we refined the or- 
bital parameters of this multi-planet system and derived a new transit ephemeris 
for HD168443b. The reduced uncertainties in the transit window make a pho- 
tometric transit search practicable. Photometric observations acquired during 
predicted transit windows were obtained on three nights. CTIO 1.0 m photome- 
try acquired on 2010 September 7 had the required precision to detect a transit 
but fell just outside of our final transit window. Nightly photometry from the 
T8 0.8 m Automated Photometric Telescope (APT) at Fairborn Observatory, 
acquired over a span of 109 nights, demonstrates that HD 168443 is constant on 
a time scale of weeks. Higher-cadence photometry on 2011 April 28 and June 25 
shows no evidence of a transit. We are able to rule out a non- grazing transit of 
HD168443b. 

Subject headings: planetary systems - techniques: photometric - techniques: 
radial velocities - stars: individual (HD 168443) 



Introduction 



The number of known exoplanets has grown rapidly in the last decade, with over 600 
confirmed exoplanets known to dat^. While most of these planets have been discovered 
using radial velocity techniques, the number of known transiting planets has increased 
signifi cantly due to dedicated transit s urveys like the space-based Kepler ( iBorucki et al. 
201 ih and CoRoT | JSarge et al. 2008h missions, and ground-based transit searches like 



t he Hungarian Autona ated Telesc ope Network (HATNet)! J Sakos et al. 2004h . Super WASP 
I ( PoUacco et al. 2006 ). and XO ( McCuUough et al. 2005). The price for efficient opera- 
tion of these wide-field transit surveys, though, is that most of the candidate stars tend to 
be fainter than those being surveyed by radial velocity. Of the over one hundred transiting 
planet host stars known, the sample of bright stars (V < 9) wi t h transiting planet s is still lim- 
i ted to only nine sta rs: HD20945 8 ffcharbonneau et al.ll200ol : iHenrv et al.lboooh HD189733 
teouchv et all2005[ ). HD149026 dSato et al.ll2005l ). HD17156 JSarbieri et al.ll2007h . HD80606 



exoplanet.eu 
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dMoutou et aPboodlFossev et al.lboOQl ). HD97658 jHenrv et al.lboill ) and 55 Cnc e| Jwinn et al. 201 ih 
(all of which were discovered by radial velocity surveys), while WASP-33b and HAT-P-2b 
were discovered by transit surveys and confirmed by radial velocity followup. 

The discovery of additional bright transiting planet hosts is advantageous in further en- 
abling studies of the atmospheric constituents of exoplanets. Even with the largest ground- 
based telescopes, transmission spectroscopy to probe the atmospheres of these exoplan- 
ets has, largely, been accomplishe d only for the bright est t argets. Using high- resolution 
spectroscopy on 8-10 m telescopes, iRedfield et al.l fcOOSl ) and [Snellen et al.l (120081 ) detected 
sodium absorption in the transmission spectra of HD189733b and HD209458b, respectively. 



More recen tly, using t h e nar rowband tunable filter imager on the 10 m GTC, I Colon et al. 



( I2OIOI ) and ISing et al.l ( 1201 ll ) have detected the signature of potassium absorption in the 
atmospheres of HD 80606b and X0-2b. 



Our ongoing Transit Ephemeris Refinement and Monitoring Survey (TERMS. iKane et al 



(I2OO9I )) project focuses on bright stars {V < 9) with known exoplanets and orbital periods 
greater than 10 days in an effort to refine the orbital parameters with additional radial ve- 
locity observations and then observe the targets photometrically within their revised transit 
windows. Transits detected around such bright stars would provide perfect candidates for 
spectroscopic follow up. In addition, with periods greater than ten days, the planet popula- 
tion searched by TERMS is not easily du plic able by ongo i ng gro und-based transit surveys as 
demonstrated by I von Braun et all (120091 ) . In iKane et al.l (I2OIII ) we presented the ephemeris 
revision and the search for a transit around HD 156846. In this paper, we present additional 
radial velocities and refine the transit ephemeris for the bright star HD 168443, which is 
known to have multiple companions. We present new photometry that allows us to rule out 
transits of HD 168443b. 



2. HD 168443 



HD 168443 (GJ 4052, HIP 89844, TYC 5681-1576-1) is a bright {V = 6.92 ) 05 dwarf 
known to possess two substellar companions, forming a dyna mically active system (IVeras fc Armitage 
2007h . HD 168443b jMarcv et al]ll999l : IWright et allbood ) has a reported Mpsini = 7.8 ± 
0.259Mjup, an orbit with a period of 58.11 days and a large eccentricity of e = 0.53. 
HD 168443c, a brown dwarf companion (lUdry et al.ll2002l : IWright et al.ll2009l ). has MpSim = 
17.5 ± 0.6 5Mj,;.p, an or b ital p eriod of ~1748 days and a moderate eccentricity e = 0.21. 
Using the Ivan LeeuwenI (120071 ) re-reduction of the Hipparcos data, iReffert fc Ouirrenbach 



(I2OIII ) derive a mass of SO.SZli^Mjup and a 3a upper mass limit of 65Mj„p, confirming 



that this object is indeed substellar. The 3a lower limit does not exclude an inclination of 
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90°, so the minimum mass derived from t he radial velocities applie s. From CORALIE ra- 
dial velocities and the Ivan LeeuwenI (|2007| ) Hipparcos re-reductions, ISahlmann et al.l (1201 if ) 
concluded that, while their formal solution for the mass of HD 168443c matched that of 
Reffert Sz Quirrenbach the mass was of "low confidence." They are unable to set an 

upper mass limit because the radial velo city orbit is n o t fully covered by their CORALIE 
obsevations. Dynamical simulations by IVeras fc FordI (120101 ) show that almost no stable 
systems can exist for mutual inclinations between HD 168443b and c of 60°-120°. 

The eccentric orbit of HD 168443b increases its transit probability significantly above 
what one would expect for a planet in a circular orbi t with the sa r ne per iod. The new orbital 
parameters, along with the formalism outlined in iKane et al.l ( 120091 ). result in a transit 
probability of 3.7% compared to 2.5% for a circular orbit. While th e atmospheric scale 
heigh ts for massive, relatively cold planets are expected to be small (jVidal-Madjar et al. 
201ll ). the brightness of HD 168443 (more than twice as bright as HD 209458, the third 
brightest star known to have a transiting planet) makes such detections possible with large 
ground-based telescopes. 

The discovery of transits of the inner planet would also constrain the possible inclinations 
of the outer brown dwarf companion, enabling additional dynamical investigations. The 
predicted transit probability of 3.7%, coupled with the fact that HD 168443 is a very bright 
star, makes it an intriguing target in our ongoing attempts to discover long-period transiting 
planets. 



3. Stellar Properties 



We used Spectroscopy Made Easy (SME: IValenti fc Piskunovlll996l ) to fit high-resolution 
Keck-HI RES spectra of HP 16844 3, applying the wavelength intervals, line data, and method- 
ology of IValenti fc Fischerl (120051 ). We further constraine d the surface gravity using Yonsei- 



Yale (Y^) stellar str ucture models (IDemarque et al 



20041) and revised Hipparcos parallaxes 



(Ivan LeeuwenI I2OO7I ). with the iterative method of IValenti et al.l (120091 ). The resulting stellar 
parameters listed in Table [1] are effective temperature, surface gravity, iron abundance, pro- 
jected rotational velocity, mass, and radius. HD 168443 li es 1.05 mag ab ove the Hipparcos 
average main sequence (My versus B — V) as defined by IWrightl (120051 ). These properties 



are consistent with an evolved metal-rich G5 star. The stellar radius, i?^ = 1.51±0.06 Rq, is 



crucial for estimating the depth and duration of a planetary transit. This value is consiste nt 
with the interferometrically measured = 1.58 ± 0.06 Rq (Ivan Belle fc von Braunll2009l ). 



In addition, we computed the level of stellar activity in HD 168443 from the strength of 



- 5 - 



the Ca II H &: K lines, which give cahbrated Shk values on the Mt. Wilson scale and log-R^K 
values (jlsaacson fc Fischer 2010 ). The median of log-R^j^ and ^hk values are listed in Table 
[T] and demonstrate that HD 168443 is chromospherically quiet. Additional examination of 
the available history of Ca II H & K measurements show no significant long-term variation 
in 5'hk- 



4. Keck-HIRES RV Measurements and Revised Orbital Model 



4.1. Measurements 



We observed HD 168443 using the stan dard procedures of the California Planet Search 
(CPS) for the HIRES echelle spectrometer (jVogt et al.lll994l ) on the 10 m Keck I telescope. 
These measurements span fifteen years, 1996 July to 2011 March, and comprise one of 
the longest RV datasets presented for a star with one or more known p lanets. The initial 



measurements in this time series were used to dis cover the two plane ts (IMarcy et al.lll999 



200ll ) while later measurements refined the orbits (jWright et al.l 120091 ). The full set of mea- 
surements presented here refines the orbit further and gives an accurate predicted transit 
ephemeris with which we search for photometric transits. 

The 130 Keck RV measurements (Table [2]) were made from observations with an iodine 
cell mounted directly in front of the spectrometer entrance slit. The dense set of molec- 
ular absorption lines imprinted on the stellar spectra provide a robust wavelength fiducial 
against which Doppler shifts are measured, as well as strong const raints on the shape of 



the spectrometer in strumental profile at the time of each observation (IMarcy fc Butlerlll992 



Valenti et al.lll995l ). We measured the Do ppler shift of each star-times-iodine spectrum us 



i ng a modeling procedure descended from iButler et al.l (119961 ) as described in iHoward et al. 



(120091 ). The times of observation (in barycentric Julian days), relative RVs, and associated 
errors (excluding jitter) are listed in Table [21 We also observed HD 168443 with the iodine 
cell removed to construct a stellar template spectrum for Doppler modelling and to derive 
stellar properties. 



4.2. Keplerian Model 



We modelled the Keck RVs as the superposition of the Keplerian interactions from two 
planets with the star, plus a linear trend in velocity due to a distant and mass ive third 
companion. We used the orbit fitting techniques described in iHoward et al.l (120101 ) and the 
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partially linearized, least-squares fitting procedure described in IWright fc HowardI (120091 ) . 
Each velocity measurement was assigned a weight, constructed from the quadrature sum 



of the measurement uncertainty (ctrv) and a jitter term (cxjitter), i-e. 



w 



jitter 



)• 



We chose jitter values of cxjitter = 3 and 2 ms~^ for measurements before and after the HIRES 
upgrade in August 2004. These values are consis tent with the expected jitter of a slightly 
evolved early G star observed with Keck/HIRES IWrightl (120051 ) . Possible sources of jitter 
include stellar pulsation, rnagnet i c cycles, gran ulation, undetected planets, and instrumental 
effects Jlsaacson fc FischeilEoiol : Iwrightlboosh . 



Our best-fit orbital model is presented in Table [3] and Figure [H The Keplerian parameter 
uncertainties for each planet were derived using a Monte Carlo method (IMarcy et al.ll2005l ) 
and account for correlations between parameter errors. Uncertainties in M sin i and a refiect 
uncertainties in M^, and the orbital parameters. We considered and rejected more compli- 
cated models having a third planet and/or a quadratic velocity trend because of statistically 
insignificant changes in compared to the adopted model (Table [3]). 



5. Transit Ephemeris Refinement 

The revised orbital solution presented in Table Ofor this multi-planet system, along with 
the stellar properties in Table [H allow us to constr uct an accurate tran s it eph emeris from 



which to conduct a search for transits. As shown by iKane fc von BraunI (120081 ). the transit 
probability of a planet is intricately related to both the orbital eccentricity and the argument 
of periastron. Figure [2] depicts the orbits of the planets relative to the observer line-of-sight 
and shows how the eccentricities and orientations affects the star-planet separation along 
that line. 



We use the models of iBodenheimer et al.l (120031 ) to estimate a radius for HD 168443b 



of Rp = l.llRjup, which takes into account orbital parameters and the stellar flux received 
by the planet. This results in a transit probability of 3.7% and a predicted transit depth 
of 0.6%. The predicted transit duration is 0.36 days, or 8 hours and 40 minutes. The 
time of mid transit shown in Table [3] is for 2011 March 1. The calculation was performed 
using a Monte-Carlo bootstrap which propagates the uncertainties in the orbital parameters 
forward to the time of transit. The uncertainty in this predicted time is small; only ~ 35 
minutes. Thus, the predicted duration of the transit window for this date is 9 hours and 
50 minutes (the sum of the predicted duration ± 1-a deviation), overwhelmingly dominated 
by the predicted transit duration rather than the uncertainty associated with the orbital 
parameters. The predicted duration makes complete coverage of the transit window from a 
single ground-based longitude difficult, though the predicted depth is sufficient to rule out 
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transits from observations during times of ingress or egress only. 



6. Hipparcos Photometry of HD 168443 



The Hipparcos missi on o bserved the brig 



Robichon fc Arenoul (|2000| ) and ICastellano et al. 



i test s tars in the sky over many epochs. 
J2000h detected the transi t of H P 209458b 



in the Hipparcos epoch photometry, and lHebrard fc Lecaveher Pes Etangd (120061 ) were able 
to detect multiple transits of HP 189733b, leading to a significant improvement in the deter- 
mination of its period. Hipparcos photometry has therefore been demonstrated to be precise 
enough to detect the transit of a Hot Jupiter around stars that are fainter than HP 168443 
(V=6.92). The top panels of Figure[3]show the Hipparcos photometry of HP 209458 plotted 
against Julian Pate and also against orbital phase of the star's Hot- Jupiter companion. A 
few observations fall within the modern transit window and show a clear dimming. The 
bottom two panels are similar plots f or HP 168443. Using our n ew radial velocity observa- 
tions, we followed the prescription of iRobichon fc Arenoul (120001 ) to create the phased plot 
with the predicted transit set at zero phase. Three Hipparcos photometric observations 
acquired around BJP 2448050 lie inside of our predicted transit window. Unlike the case 
for HP 209458, the Hipparcos photometry show no evidence for a transit. However, the ex- 
pected transit depth for HP 168443b is significantly smaller than HP 209458; the Hipparcos 
photometry alone would have been an unreliable guide to prove or disprove the occurrence 
of transits in HP 168443. Nevertheless, we advocate such a check for bright stars with 
sufficiently precise orbital ephemerides. 



7. Photometry of the Revised Transit Window 

In addition to the new TERMS photometry presented below, we also observed HP 168443 
with the Southeastern Association for Research in Astronomy (SARA) 0.6 m telescope at 
Cerro Tololo Inter-American Observatory. However, due to poor weather conditions, the 
SARA-S measurements exhibited scatter that was significantly higher than the predicted 
transit depth. Therefore, we do not discuss these data further, but we note that planning 
such multi-site observations is necessary in our TERMS search for long-period transits. 
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7.1. Cerro Tololo Inter- American Observatory (CTIO) 

Before we had access to the lates t Keck-HIRES radial velocities, we used published 



orbital parameters ( iWright et al.ll2009l ) to calculate the transit ephemeris and to schedule 
an observing run for 2010 September 7with the CTIO 1.0 m telescope and Y4KCam CCD 
detector. The observations were made through a Johnson- Morgan R-band filter; instru- 
mental magnitudes of HD 168443 and com parison stars were extracted from the images 



with an IDL implementation of DAOPHOT (jStetson P. 1119871 ). Relative fluxes (Everett & 



Howell 2001) of HD 168443 were computed with respect to the two stable comparison stars 
(TYC 5681-1450-1 and TYC 5681-1458-1). 

The results are plotted in Figure 4. The solid curve represents the predicted transit 
fluxes computed from our new orbital elements in Table 2, assuming on-time, central transits 
with a predicted depth of 0.6% or 0.006 flux units. The scatter of the measurements is 0.0027, 
easily sufficient to detect the predicted transits. However, these measurements cover only 
the later part of the predicted transit window; an egress late by la and 3 a are represented 
with dashed blue and red lines, respectively. The CTIO data show only that late transits do 
not occur. 

We have presented these measurements because CTIO photometric monitoring is an 
essential component of the TERMS strategy. Had the ephemerides (based on orbital param- 
eters in the literature) been more precise, the CTIO photometry has the requisite precision 
to detect a predicted on-time transit. Our results demonstrate the precision achievable from 
CTIO with a typical TERMS target. 



7.2. Fairborn Observatory 

We obtained additional photometric observations of HD 168443 with the T8 0.8 m au- 
tomatic photometric telescope (APT) at Fairborn Observatory in southern Arizona. The 
T8 APT uses a two-channel precision photometer with two EMI 9124QB bi-alkali photo- 
multiplier tubes to make simultaneous measurements in the Stromgren b and y passbands. 
The telescope was programmed to make nightly differential brightness measurements of 
HD 168443 with respect to the comparison star HD 166664. Three consecutive differential 
measurements were co-added to create a single nightly differential magnitude. To improve the 
precision of these brightness measurements, we combined the individual b and y differential 
magnitudes into a mean {b + y)/2 "pass band" . The typical precision of a single observation 



on good nights is ~ 0.0015 mag. See (lHenrylll999l ) for further details on telescope design 



and operations, data reduction, cahbrations, and data precision. 
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Between 2011 March 2 and June 19, the APT collected 107 nightly observations of 
HD 168443 with respect to HD 166664; the differential {b + y)/2 magnitudes were converted 
to relative fluxes and are plotted in Figure 5. The nightly observations scatter about their 
mean flux, indicated by the dashed line in Figure 5, with a standard deviation of 0.0016, 
which is consistent with constant stars. Periodogram analyses from one to 100 days reveal 
no significant periodicity. We conclude that HD 168443 is constant on its rotation timescale. 
A least-squares sine fit of the nightly observations on the orbital period of 168443b gives a 
semi-amplitude of 0.00026 ± 0.00023 flux units. This very low limit to brightness variability 
in HD 168443 on the 58-day orbital period indicates that r otational modulati on of starspots 



is not the cause of the radial velocity variations, see e.g., ( IQueloz et al.ll200l[ ) 



We also used the T8 APT to monitor HD 168443 on two nights in the 2011 observ- 
ing season when our new ephemeris predicted additional transit events. Each monitoring 
observation consists of a single differential measurement rather than the mean of three ob- 
servations, as we used for the nightly observations; thus the monitoring observations will 
have more scatter than the nightly observations. 

On 2011 April 28 UT, an observable egress was to occur at BJD 2,455,679.937. We 
successfully monitored the star and obtained the 80 measurements plotted in Figure 6. The 
star was still two months before its opposition, so the start time was delayed until the star 
rose above an airmass of ~ 2.0; the observations ended 3.3 hours later (at dawn). The 
sudden increase of the scatter after predicted egress is the result of a plume of smoke from 
one of the many wildfires burning throughout southern Arizona at the time. The standard 
deviation of the entire data set is 0.0040 flux units. The solid curve in Figure 6 represents 
the predicted fluxes for an on-time, central transit. The dotted blue lines represent the ± la 
uncertainty in the predicted time of central transit. The difference between the mean flux of 
the 47 pre-egress observations and the 33 post-egress observations is only 0.0005 ± 0.0011. 
Given the tight limit on the transit time (35 min or 0.024 day), these observations rule out 
central transits with the predicted depth of 0.006 with a SNR of ~ 5:1. 

The next predicted transit was calculated to occur 58 days later on 2011 June 25 UT, 
centered at BJD 2,455,737.862. HD 168443 was at opposition during this transit, so we were 
able to acquire 165 observations over an interval of 7.1 hours. The declination of HD 168443 
is approximately —10°, so the airmass values at the start and at the end of the night were 2.25 
and 2.51, respectively, for the the east and west observing limits. Atmospheric extinction 
was also significantly higher than normal at the beginning of the night due to airborne dust, 
which gradually settled out over the course of the night. Therefore, in addition to reducing 
the transit observations with larger than normal extinction coefficients, we also removed a 
linear trend of approximately 0.008 flux units via a linear least-squares fit. The residuals 
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from the line fit are plotted in Figure 7, again compared with predicted transit fluxes. The 
overall scatter of the 165 observations is 0.0039 flux units. 

Finally, we estimate an upper limit to the possible transits of HD 168443b. The mean 
differential magnitudes of the T8 nightly observations in Figure 5 and of the mean of the first 
and second T8 transit monitoring observations in Figures 6 & 7, before they were converted 
to relative fluxes, are -0.24041 ± 0.00017 mag, -0.23961 ± 0.00050 mag, and -0.24044 ± 
0.00037 mag, respectively. So, the observed difference between the mean of the nightly 
observations and the mean of the transit observations is only 0.00038 ± 0.00045 mag or 
0.00035 ± 0.00041 flux units. Therefore, we should have been able to detect a transit depth 
of 0.2% or 0.002 at a SNR of almost 5:1, but our observations show no evidence of a transit. 

7.3. Full Coverage of the Transit Window 

The entire transit window was covered during the three monitoring nights. Initially, we 
were able to rule out only a late egress using 1-meter CTIO photometry from 2010 September 
8 with (3-cr) confldence. We further conflrmed this result with the APT photometry on June 
25 by ruling out a late ingress (3-0"). The early (l-a) and on time egress was ruled out with 
APT photometry from April 28, 201. In Figure Owe present a phase plot containing all thee 
monitoring nights. Since both data sets from the T8 APT used the same reference star, we 
placed the median of the measurements on the same y-axis scale. The CTIO data set was 
adjusted with an offset using the overlapping points from CTIO and the APT. The center 
of predicted transit is at phase 0.0, marked by a solid vertical line. The dotted blue lines 
represent 1-a early and late windows while the dashed red lines represent 3-0" deviation from 
the center of the predicted transit. 

8. Discussion 

As part of our ongoing Transit Ephemeris Reflnement and Monitoring Survey (TERMS), 
we present revised orbital parameters for the HD 168443 system, based on 130 radial velocity 
measurements with Keck-HIRES that span almost 15 years. Using the transit ephemerides 
derived from the revised orbital parameters, we searched for a transit using telescopes from 
CTIO, Fairborn Observatory, and SARA-S. We flnd no evidence of a detectable transit. The 
presence of a non-grazing transit corresponding to our model (l.li?j planet) is ruled out 
at a high level of confldence with high precision photometry acquired by the APT. Grazing 
transits or transits with a planet radii as small as 0.58 Rj^p (which would yield densities 
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much too high) are formally excluded at 1-cr confidence, though the smaller duration of 
such transits implies that the time of our photometric observations could have missed the 
ingress and egress. Using our orbital solution, and the planetary and stellar radii presented 
in this paper , we derive an upper limit on inclination of the system at 87.8 ° using methods 
described by iKane fc von BraunI (120081 ). 



Even with a number of RV observations, determining precise orbital parameters for one 
component of a multi-planet sys tem can be difficult. The recent di scovery of transits of 
55 Cn c e by | ( Winn et al. 2011 ). based on a new orbital period by Dawson fc Fabry cky 



(I2OIOI ). illustrates the insidious impact of aliases and harmonics. Nevertheless, these effects 
can be mitigated by careful observation and analysis, and the value of additional precision 
radial velocities and revised orbital parameters cannot be overstated. 

While we do not see transits in the HD 168443 system, the experimental approach 
outlined here, combining high precision radial velocity with multiple photometric telescope 
facilities, is worth pursuing in the quest for new transiting planets around bright nearby 
stars. 
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Table 1. Stellar Properties 


Parameter 


Value 


Reference 


Mv 


4.198 


van Leeuwen {2QQ^^ 


B-V 


0.724 


Bessell ('20001 


V 


6.92 


Bessell (2000) 


Distance (pc) 


37.4 ± 1.0 


van Leeuwen f2007) 


Teff (K) 


5491 ± 44 


This work 


logs 


4.07 ±0.06 


This work 


[Fe/H] 


+0.04 ±0.03 


This work 


f sini (kms^"'^) 


2.20 ±0.50 


This work 


(Mo) 


0.995 ±0.019 


This work 


R* [Rq) 


1.51 ±0.06 


This work 


l0gi?iiK 


-5.088 


This work 


S'hk 


0.148 


This work 
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Table 2. Keck Radial Velocities 





Radial Velocity 


Uncertainty 


BJD - 2440000 


(ms-i) 


(ms-i) 


10276.90890 


-309.29 


1.73 


10603.01184 


-33.47 


1.02 


10665.86781 


-69.33 


1.24 


10713.73770 


-73.02 


1.14 


10714.76649 


-71.38 


1.13 


10955.01039 


-5.57 


1.12 


10955.95862 


-5.13 


1.02 


10957.07105 


-4.68 


1.05 


10981.88012 


-568.69 


1.01 


10982.89132 


-492.38 


1.10 


10983.07690 


-472.56 


1.10 


10983.82231 


-420.40 


1.16 


10984.06138 


-407.63 


1.14 


11009.87009 


22.04 


1.34 


11010.05994 


24.94 


0.93 


11010.85123 


20.67 


1.17 


11011.86077 


25.19 


1.37 


11012.95413 


19.19 


1.14 


11013.06816 


20.31 


0.82 


11013.82791 


15.69 


1.15 


11013.92983 


15.33 


1.20 


11042.95557 


-326.88 


1.09 


11043.95602 


-276.32 


1.19 


11050.81406 


-77.68 


1.26 


11068.77042 


54.21 


1.08 


11069.78596 


60.86 


1.09 


11070.79807 


52.81 


1.09 


11071.76998 


50.98 


1.11 


11072.76271 


47.77 


1.26 


11074.78514 


39.07 


1.09 


11228.16111 


91.81 


1.07 


11229.14942 


95.62 


1.32 


11311.04174 


151.35 


1.27 


11312.07757 


137.66 


1.23 


11313.07589 


124.69 


1.18 


11314.08699 


103.43 


1.20 



Table 2 — Continued 



Radial Velocity Unccirtainty 
BJD- 2440000 (ms"!) (ms-i) 



11341.02661 


109.66 


1.11 


11341.90588 


129.71 


1.03 


11342.97144 


142.29 


1.20 


11367.86065 


214.31 


1.30 


11368.84314 


205.20 


1.40 


11370.01150 


188.83 


1.14 


11370.91933 


171.14 


0.73 


11371.91024 


153.12 


1.11 


11372.87811 


131.18 


1.16 


11373.79496 


105.31 


1.44 


11409.85005 


276.28 


1.14 


11410.84584 


278.53 


1.29 


11411.84641 


288.70 


1.21 


11438.73808 


-282.82 


1.19 


11439.73057 


-414.58 


1.18 


11440.71830 


-527.69 


1.22 


11441.73997 


-616.71 


1.20 


11679.04801 


-109.80 


1.09 


11680.07181 


-20.27 


1.29 


11703.01759 


488.68 


1.19 


1 1 703 98840 


488 QO 


1 17 


11705.03899 


493.45 


1.27 


11705.95846 


499.46 


1.30 


11707.08125 


502.77 


1.20 


11754.85883 


473.14 


1.66 


11755.90986 


489.51 


1.06 


11792.74760 


-338.38 


1.01 


11793.80008 


-223.10 


1.25 


11882.68333 


496.02 


0.88 


11883.68272 


498.30 


0.86 


11983.15885 


333.08 


1.01 


11984.15396 


335.98 


1.35 


12004.12317 


357.01 


1.33 


12005.14737 


352.21 


1.37 


12007.13381 


329.80 


1.16 


12008.03597 


311.14 


1.37 



Table 2 — Continued 



Radial Velocity Unccirtainty 
BJD- 2440000 (ms"!) (ms-i) 



12009.10491 


294.16 


1.24 


12030.98196 


-27.97 


1.26 


12061.94969 


278.32 


1.37 


12062.96034 


259.31 


1.54 


12094.88634 


78.37 


1.40 


12096.93360 


110.11 


1.52 


12098.01419 


126.96 


1.38 


12099.02269 


145.64 


1.34 


12099.93319 


147.47 


1.23 


12100.94388 


158.52 


1.41 


12101.88510 


171.72 


1.42 


12127.86347 


57.89 


1.51 


12128.79403 


30.56 


1.23 


12133.78671 


-210.61 


1.53 


12160.80699 


79.11 


1.48 


12189.76875 


-156.78 


1.43 


12445.93556 


-192.60 


1.38 


12486.81165 


-930.27 


1.38 


12515.75755 


-85.29 


1.44 


12536.74799 


-258.80 


1.39 


12572 69122 


-74.58 


1.31 


12713.14689 


-280.58 


1.33 


12778.02605 


-891.12 


1.36 


12804.06082 


12.02 


1.25 


12834.88152 


-730.70 


1.31 


12848.80396 


-133.69 


1.29 


12855.96590 


0.22 


1.28 


12898.71740 


-579.54 


1.16 


13154.04176 


260.17 


1.49 


13180.90222 


-183.23 


1.37 


13195.84965 


58.60 


1.13 


13238.88150 


-123.80 


1.17 


13301.73743 


-565.16 


1.16 


13546.89064 


334.18 


1.14 


13842.12630 


125.12 


1.18 


13927.87788 


-3.76 


1.15 
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Table 2 — Continued 





Radial Velocity 


Uncertainty 


BJD - 2440000 


(ms~^) 


(ms~^) 


13984.83683 


-51.21 


0.96 


14314.99539 


-101.88 


1.26 


14335.95562 


-191.79 


1.23 


14343.88345 


-620.98 


1.00 


14344.94159 


-744.93 


1.11 


14398.74642 


-344.38 


1.08 


14546.11545 


-16.53 


1.17 


14548.15404 


10.39 


1.27 


14720.84072 


102.86 


1.11 


14956.12833 


289.18 


1.53 


14985.11014 


-473.97 


1.33 


15014.97707 


332.37 


1.13 


15026.96186 


318.51 


1.22 


15106.74809 


-279.16 


1.20 


15286.11499 


173.04 


1.21 


15322.08481 


403.90 


1.21 


15343.00253 


103.85 


1.36 


15374.84199 


453.14 


1.16 


15378.82083 


413.14 


1.21 


15403.81739 


198.38 


1.22 


15490.73506 


346.47 


1.30 


15636.13540 


-128.87 


1.05 
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Table 3. Keplerian Orbital Model 



Parameter 



Value 



HD 168443 b 
P (days) 

Tc" (JD - 2,440,000) 
Tp^ (JD - 2,440,000) 

e 

K (ms-^) 
u (deg) 
Msim (Mjup) 
a (AU) 

HD 168443 c 
P (days) 

Tc" (JD - 2,440,000) 
Tp^ (JD - 2,440,000) 
e 

K (ms-i) 
w (deg) 
Msini (Mjup) 
a (AU) 

System Properties 
7 (ms~-^) 

di;/di (ms~''^yr~^) 
Measurements and Model 



obs 



rms (ms ) 



58.11247 ±0.0003 
15621.637 ±0.0156 
15626.199 ±0.024 
0.52883 ±0.00103 
475.133 ±0.9102 
172.923 ±0.139 
7.659 ± 0.0975 
0.2931 ± 0.00181 



1749.83 ±0.57 
15599.9 ± 1.187 

15521.3 ±2.2 
0.2113 ±0.00171 
297.70 ±0.618 
64.87 ±0.5113 

17.193 ±0.21 
2.8373 ± 0.018 



-46.533 ± 0.552 
-0.00868 ± 0.00025 



130 
3.90 
1.44 



^Time of transit. 

^Time of periastron passage. 
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Fig. 1. — RV measurements of HD 168443 from Keck-HIRES (filled circles) with Keplerian 
orbital model (dashed lines). The top panel shows the RVs phased to the orbital period of 
HD 168443 b with the model for the other planet and linear trend subtracted. Open circles 
represent the same RV measurements wrapped one orbital phase. The bottom panel shows 
the RV time series illustrating the variations due to HD 168443 c, with the linear velocity 
trend and the orbit of HD 168443 b subtracted. 
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V to observer 



Fig. 2. — The orbits of the planets HD 168443b and HD 168443c shown in sohd hnes. Orbits 
of Mercury (eccentricity set to zero for clarity) through Jupiter are represented with dashed 
ovals. 
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Fig. 3. — Hipparcos Photometry. Top left: Photometric measurements of HD 209458 over 
the course of three years. Top right: the same data plotted vs. orbital phase, with the 
transit corresponding to zero phase. The transit is clearly visible and is marked by two 
vertical lines. Bottom left: Hipparcos Photometric measurements of HD 168443 taken over 
the course of three years. Bott o m rig ht: Similar to the plot for HD 209458 (following the 
method of iRobichon fc Arenoul (120001 )). The predicted window is bounded by two vertical 
lines, with three Hipparcos measurements obtained during the predicted transit. No transit 
is evident. 
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Fig. 4. — Relative photometry from CTIO. The standard deviation of 0.0027 in the nor- 
mahsed flux is precise enough to detect the predicted transit, which would causes a 0.006 
decrease. The solid line represents the predicted transit after ephemeris refinement using 
Keck HIReS data, with the predicted egress occurring right before the measurements were 
acquired. The dotted blue and the dashed red lines represent the 1-cr and 3-a uncertainties 
in the time of egress. 
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Fig. 5. — Photometry from the T8 APT consisting of 107 measurements over the span of 
109 days. The dashed line represents the normalised flux. The standard deviation from the 
mean is 0.0016. 
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Fig. 6. — Photometric observations of HD 168443 acquired with the T8 APT at Fairborn 
Observatory during the predicted transit of 2011 April 28. The standard deviation of 0.0040 
flux units is sufficient to detect the predicted transit egress (sohd fine) if present. The dotted 
blue lines represent the ± 1 — cr uncertainty in the transit time. There is no evidence for an 
egress event. 
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Fig. 7. — Photometric observations of HD 168443 with the T8 APT during the predicted 
transit of 2011 June 25. The standard deviation of 0.0039 flux units is sufficient to detect the 
predicted transit egress (solid line) if present. The dotted blue lines represent the ± 1 — cr 
uncertainty in the transit time. Again, there is no indication of a transit. 
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Fig. 8. — Phase diagram of the observations from all three nights of transit monitoring. The 
filled and unfilled circles are the T8 APT measurements from 2011 June 25 and April 28, 
respectively. The filled triangles are the CTIO observations from 2010 September 7. As 
in earlier figures, the solid line represents the predicted flux changes during a transit. The 
dotted blue and the dashed red lines represent ±1 — 0" and ± 3 — cr deviations in the time 
of transit, respectively. 



